Functional renormalization group for quantized anharmonic oscillator 
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Functional renormalization group methods formulated in the real-time formalism are applied to 
the O(N) symmetric quantum anharmonic oscillator, considered as a + 1 dimensional quantum 
field-theoric model, in the next-to-leading order of the gradient expansion of the one- and two- 
particle irreducible effective action. The infrared scaling laws and the sensitivity-matrix analysis 
show the existence of only a single, symmetric phase. The field-independent term of the wavefunction 
renormalization turned out to be negligible, but its field-dependent piece is noticeable. It is shown 
that the infrared limits of the running couplings depend on the renormalization group scheme used, 
when the perturbation expansion in the bare quartic coupling is truncated keeping the terms up to 
CD ■ the second order. 
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I , Our main goal in this paper is to investigate the role of wave-function renormalization in the quantum properties of 
£S| ' the quartic anharmonic oscillator considering it as a O(space)+l(time)-dimensional quantum field-theoric model, and 
applying various internal space renormalization group (RG) schemes to the one- (1PI) and two-particle irreducible 
(2PI) effective actions in the truncated gradient expansion which goes beyond the local potential approximation (LPA) 
with the incorporation of wave-function renormalization terms. In that sense it is a continuation of the work presented 
in [1]. On the one hand, we apply a Callan-Symanzik (CS) type internal space (IS) RG scheme to the 1PI effective 
' action, where an imaginary mass parameter serves as the control parameter of the RG evolution. On the other hand, 
another IS RG scheme is applied to the 2PI effective action, where the bare quartic coupling gs represents the control 
l— ~~ '■ parameter. In both cases, the gradient expansion is used keeping the wavefunction renormalization with, and the 
I | local potential is truncated beyond its quartic term. The results obtained in the framework of various RG schemes in 
■ the second-order perturbation expansion in the coupling gs are also compared. 
7-H ' The quantized anharmonic oscillator is an elementary toy model representing a building block of more sophisticated 
and/or realistic models finding several physical applications. Let us mention a few of them without pretending to 
completeness. The investigation of quantum properties of nanomechanical oscillators has became of interest in order 
to understand the quantum behaviour of nanodevices [2] . The statistical physical properties of chains of anharmonic 
Q\ • oscillators, that of a one-dimensional Ginzburg-Landau system with quartic anharmonicities is of great interest in solid 
CD state physics. 2- and 3-dimensional aggregates of such chains describe anisotropic anharmonic solids (e.g. ultrathin 
ferroelectric films). Putting such systems in inhomogeneous external field enables one to manipulate individual 
monomers, making such films very attractive for molecular electronics [3]. 1-dimensional Ginzburg -Landau systems 
. are very often investigated by means of the transfer matrix method in which the evaluation of the free energy of an 
N + 1 dimensional classical system is equivalent to the evaluation of the ground state energy of an TV-dimensional 
quantum system, the so-called dual counterpart of the classical system. In general this dual problem is solved to 
understand its classical counterpart [4] , [5] . Another field where the anharmonic oscillator model has got application 
is the study of the quantum physics of systems showing up chaotic behaviour classically, ie. the study of quantum 
chaos [6], [7], [8]. Also the possibility of quantum effects induced chaos has recently been debated on the problem of 
the quantized Duffing oscillator [9]. 

These widespreading possibilities of its applications explain the efforts to understand the quantum physics of the 
anharmonic oscillator in all detail and to invent various theoretical tools for its investigation. Already in the pioneering 
works [10], [11] it has been pointed out that the Rayleigh-Schrodmger perturbation series for the ground state energy 
of the quartic anharmonic oscillator diverges. Since then much efforts have been made to overcome that problem 
and determine the eigenvalues and eigenfunctions of the anharmonic oscillator by various methods like the strong- 
coupling expansion [12], the integration of operator differential equations [13], multi-scale perturbation theory [14], 
variational methods [15], iterative Bogoliubov transformations [16], the eigenvalue moment method [17], the improved 
Hill-determinant method [18], the optimized perturbation expansion [19], a particular iterative method based on the 
generalized Bloch-equation [20], the auxiliary field method combined with loop expansion [21], as well as the quantum 
computational method [22]. A class of non-perturbative approaches is based on various RG methods [23-25] including 
functional RG methods [26, 27], [1] which treat the quantum mechanical problem as a + 1 dimensional quantum 
field-theoric one. 

The RG method has been basically developed for the purposes of describing phase transitions and critical behaviour 
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in statistical physics and very soon extended to investigate quantum field-theoric models. The RG strategy has been 
aimed to monitor the Green functions of the elementary held variable when the quantum fluctuation modes are turned 
on gradually [28] . While in the original version of the RG method a finite fraction of modes are turned on in each 
blocking steps and one generally relies on the perturbation expansion in bookkeeping the change of the dynamics, the 
functional RG methods [29] -[32] turn on an infinitesimal fraction of modes in the subsequent blocking steps. Such an 
approach is accompanied with the usage of infinitely many vertices. The manipulation of a high number of couplings 
is carried out through their generating functional for which one obtains an exact RG equation looking like a one-loop 
equation. In order to solve that one has to project it to a suitably chosen functional subspace, e.g. by means of 
the truncated gradient expansion, but that seems to be a physically better motivated approximation scheme than 
the truncation of the perturbation expansion. The functional RG strategy has been developed for the 1PI effective 
action [31], [32] with gliding momentum cut-off. According to the original RG idea one organized the quantum 
fluctuations according to their wavelengths or momentum, ie. according to their properties in the external space. 
A natural generalization came when such an ordering has achieved according to their internal-space characteristics, 
the largeness of their amplitudes. That step has been done by developing the functional generalization of the CS 
equation [33] where the mass parameter is evolved [34] and the modes of the quantum fluctuations are turned on in 
the order of their increasing amplitudes. The CS RG scheme and the functional RG schemes with gliding momentum 
cut-off give similar scaling laws in the leading order in the UV scaling regime where the cutoff is the only scale 
parameter, but approaching the IR scaling regime the cut-off dependence and the mass-dependence may differ due 
to the occurring of additional physical scales. A further generalization of the CS RG strategy is called IS RG [35], 
when any of the parameters of the model can be used as control parameter and evolving it gradually one arrives at a 
functional RG evolution equation for the 1PI effective action. One should mention that all IS RG schemes including 
the CS RG compare various theories belonging to the gradually altering values of the control parameter, while the 
RG schemes with gliding momentum cut-off provide true scaling laws for a given theory. So long the suppression 
term generating the IS RG evolution is quadratic, the evolution equation contains one-loop integrals. But it can be 
rather complicated when the suppression term is not quadratic in the elementary field variable since the excitations 
described by the higher-loop integrals appearing in the evolution equation have rather complicated structure in single- 
particle terms. That problem can be overcome applying the internal space RG to the generating functional of some 
composite operator of the elementary field, in which the suppression term is again quadratic. In particular, in the 
case of a suppression term being quartic in the elementary field variable, one may apply the RG method to the 2PI 
effective action [35], the generating functional of 2PI 2-particle correlation functions [36, 37]. That method has the 
advantage that all expressions are formally ultraviolet (UV) finite until the generally UV divergent loop integrals are 
made explicit making an Ansatz for the parametrization of the one-particle Green function in whose terms those are 
expressed. 

As mentioned above, various functional RG methods have been recently used and proved to be powerful nonper- 
turbativc tools in treating the problem of the quantum anharmonic oscillator. In [26] the effective average action 
RG scheme [31] and in [27] the Wegner-Houghton (WH) RG method [29] have been applied to the one-dimensional 
quantum oscillator with quartic anharmonicity in order to determine the energy of the ground state and that of the 
gap between the first excited state and the ground state and both have been found to increase strictly monotonically 
with increasing bare quartic coupling gs ■ It has also been thoroughly discussed that even for a symmetric double- well 
bare potential there exists a unique ground state, so that the quantum anharmonic oscillator exhibits only a single 
phase. In the WH RG approach one integrates out the high-frequency modes of quantum fluctuations in infinitesimal 
steps above the gliding sharp momentum cut-off k. The method has the disadvantage that it disables one to go 
beyond the lowest order of the gradient expansion, the so-called LPA which has also been used in [27]. Although the 
effective average action method allows for taking into account wave-function renormalization terms in the gradient 
expansion, it has not been aimed in [26]. In [1] the functional RG method has been applied to the 1PI effective 
action for single-particle quantum mechanics using the sharp gliding momentum cut-off as control parameter of the 
evolution. Truncating the expansion of the 1PI effective action at quartic terms, RG flow equations for the 1PI two- 
and four-point vertex functions have been derived. As an example of applications, the energy of the ground state and 
that of the first excited state have been determined numerically. The momentum-dependence of the proper self-energy 
has been taken with, but the role of wave-function renormalization has not been analysed directly. The authors found 
in accordance with the finding in [27] that the Lehman-expansion of the propagator is dominated by the single pole 
corresponding to the first excited state that hints on the negligible role of the wavefunction renormalization. Our pur- 
pose is to analyse the role of wavefunction renormalization more thoroughly in terms of two different IS RG schemes, 
the CS type one (CSi) with the imaginary mass as control parameter and another one (ISg) with the quartic coupling 
as control parameter. 

First, we shall apply the CSi RG scheme to the 0(A^)-symmetric generalization of the one-dimensional quartic 
oscillator considering the oscillator coordinate q = (qi, . . . , gjv) an O(N) vector in the internal space. Due to O(N) 

symmetry the anharmonic potential can only depend on powers of the O(N) scalar q 2 = J2a=i = <7a?a- O(N) 



3 



symmetry couples any one-dimensional oscillator q a (a = 1,2,... ,7V) of the system to all of the others in a rather 
particular manner even if we restrict ourselves to a quartic potential, as we shall do throughout this work. We shall 
derive the functional CSi RG equation for the 1PI effective action, T[q) being the generating functional of the 1PI 
Green functions of the elementary field variable q. In the right-hand side of the evolution equation the inverse matrix 
under the trace shall be Neumann-expanded in its off-diagonal piece and the terms up to the quadratic ones in the 
latter only kept. The solution of that truncated evolution equation shall be looked for in the next-to- leading order of 
the gradient expansion, including field-dependent wavefunction renormalization besides the local potential. In the case 
of quartic self-interaction the field-independent wavefunction renormalization emerges as a two-loop effect, therefore 
the way to generate it by an essentially one-loop evolution equation opens via including field-dependent wavefunction 
renormalization terms, as well. We shall take with only the quartic ones of those. The RG evolution of the couplings 
and their IR values are then numerically investigated in the function of the dimension N of the oscillator. One should 
notice that the truncations used may lead to erroneous results for sufficiently large dimension N and bare coupling 
qb [38]. We shall conclude that the O(N) symmetric anharmonic oscillator with any dimension N exhibits a single 
phase and the IR couplings including those of the wavefunction renormalization depend smoothly of the dimension 
TV. 

At second, we apply the ISg RG scheme to the second Legendre-transform of the generating functional of connected 
Green- functions of the bi-local composite operator qtqt', to the 2PI effective action T[G] that is the functional of 
the single-particle propagator G. According to our results obtained in the CSi RG scheme no qualitative changes 
occur with increasing dimension N in the trivial phase structure and the IR scaling laws of the anharmonic oscillator. 
Therefore we shall restrict our further discussion to the simple anharmonic oscillator with a single degree of freedom 
in the internal space, ie. with N — 1. The 2PI effective action T[G] takes its minimum at the propagator for the 
ground state and its value at that minimum determines the energy of the ground state. Parameters of the wave- 
function renormalization and the excitation energy of the first excited state can be read of from the propagator of the 
ground state. The renormalized coupling of the quartic term of the potential is determined via the second functional 
derivative of the effective action T[G] at the ground state, the two-particle propagator of the vacuum in quantum 
field-theoric terms. The RG evolution equation for the 2PI effective action shall be derived and solved with an Ansatz 
keeping the terms up to the quadratic ones in the 2-particle propagator, called expansion in the 2-particle channel 
(E2PC). Finally we shall compare the results obtained in second order of the bare coupling g B by means of the CSi 
and the ISg RG schemes, as well as with other results taken from the literature. The dependence of the contributions 
of the order O {g B ) on the RG schemes used shall be demonstrated. 

The structure of the paper is as follows. Sect. II contains the derivation of the CSi RG evolution equation for the 
1PI effective action for the 0(N) symmetric anharmonic oscillator, the derivation of the coupled set of RG equations 
for the running couplings present in the truncated gradient expansion of the effective action, and the comparison of 
the numerical solutions obtained in the independent mode approximation (IMA), the LPA, and the approximation 
including wavefunction renormalization (AWF). The existence of only a single phase is stressed by discussing the 
sensitivity of the IR couplings to their bare values, as well. Sect. Ill starts with the derivation of the RG evolution 
equation for the 2PI effective action for the one-component (N = 1) quartic anharmonic oscillator, that is followed 
by the solutions first in the IMA, and afterwards in the truncated E2PC. Finally, the results for the various couplings 
are determined keeping the terms up to the order O (g%) and compared with the results obtained in the CSi RG 
scheme in second order perturbation expansion, as described in Appendix A3, with the WH RG scheme in the same 
approximation shortly outlined in Appendix B and the well-known results of the Rayleigh-Schrodinger perturbation 
expansion (RS PE). In Sect. IV the conclusions are drawn. The details of the derivation of CSi RG evolution equations 
are given in Appendix A 1, the loop integrals needed in both schemes are given in Appendix A 2. Appendix B presents 
a short derivation of the WH RG equation for a one-dimensional quantum mechanical system and its application to 
the anharmonic oscillator in the second order of the perturbation expansion. 

II. CSI RG FOR THE N-DIMENSIONAL OSCILLATOR 
A. RG equations 

In the Euclidean formulation of the quantum field-theoric models the CS RG scheme is realized by evolving the mass 
as a control parameter from above the UV momentum cut-off A towards its physical value continuously. The evolution 
accounts for the quantum fluctuations gradually according to their increasing amplitude. Here we realize the similar 
RG procedure in the real-time formulation by introducing an evolving imaginary mass as control parameter which 
vanishes in the IR limit. Euclidean flow equations being easier to manage numerically are then obtained by analytic 
continuation of the control parameter. The advantage of the real-time approach is that the path integral remains 
well-defined even for a quantum system with double-well potential, ie. for negative values of the bare parameter 



4 



u>q. Nevertheless one does not expect any phase transition with analytic continuation of the uJq parameter to the 
negative real axis, since there emerges no spontaneous symmetry breaking in quantum mechanics. The large amplitude 
quantum fluctuations fill up the potential wells and make it convex in the IR limit (c.f. with the findings in Refs. 
[26, 27]). 

The basic idea is that we introduce an imaginary part \i^i 2 J t q 2 into the quadratic term of the bare action Sb[q\, 
in order to suppress the quantum fluctuations in the path integral for [i 2 m A 2 . The real control parameter /j 2 is then 
continuously decreased to zero. The imaginary quadratic part ensures the convergence of the path integral during the 
whole evolution. The generator functional of the connected Green- functions of the elementary 'field variable' q t , the 
time-dependent coordinate of the linear O(N) symmetric anharmonic oscillator is given as 

Ssig] = l(lg 2 \H ^ V - ff(<Z 2 ) 2 ) ■ (1) 

For n 2 = p? B k> A 2 , the quantum fluctuations are frozen and they can be gradually 'melted out' when the control 
parameter is decreased towards /i 2 — > 0. For the sake of simplicity the mass m of the oscillator is set to unity. 
Throughout the paper we shall use the notations and conventions J t — J T ^ 2 dt, (T — > oo), J u = J^\^, (A — > 

°°)> J2a=l It fa,t9a,t = / ' 9, fa,u = J t e ZOjt f a ,t, I(a,t),(b,t') = &a,b$t,t' = S a . b S(t-t'), I( a ,w),(b,u') = <5a,6<L+w',0 = 

6a,bJ t ^e iuJt+iuj,t '6t,t>, <L=o,o = T, tiA = J2a=i St,t' S t>t' A (<*,t), («>*') = T T,a=i L A (a, u ),(a,- U ), and identical Latin 
indices denote summation like f a g a — J2a=i fa9a- The generating functional of the connected Green- functions of the 
elementary variable q is defined via the path integral 



e i w bJ = / x>qe^ SB ^+H-l (2) 



in the presence of the time-local external source j . The generating functional of the 1PI Green functions, the 1PI 
effective action as the functional of the ground-state expectation value of the coordinate, q a ^ = 5W\j]/8j a ,t is defined 
by the Legendre-transform 

-T[q] = -W\£}+i-q (3) 

for which j a .t = —5T[q]/8q a ,t, and I^ 2 ) • — —I where and r' 2 ) stand for the second functional derivatives 

of the functionals W[j) and r[q], respectively. In order to find the CS-type functional evolution equation, we evaluate 
the partial derivative of the effective action with respect to the control parameter fi 2 , 

d^T[q] = +d , 2 W\j] = +e-M \vq /UlWl^^iWH^ f S " Jw® 

~ J ~Jt 2 ~ 2 Jt OjrJa.tOjrJa^t 

2 Jt V % &3a,t$3a,t Sja,t Sja.t J 2\ H J 

Now we introduce the reduced effective action T[q] with the relation T[q] = f [q] + \i^ 2 q ■ q in terms of which the 
evolution equation takes the following form, 

d^t[q] = -^(fWy+i/i 2 )" 1 . (5) 

We shall look for the solution of the evolution equation in the truncated gradient expansion by making use of the 
Ansatz 



r M = -Tin 2 + \ I ( x V 2Uj2 ~ ^ 2 )<L<L W ~ / V( W l' W 2)^ 1 +...+ W4 ,0?o, W i5a,u- 2 g6, W3 96,W4 

^ J UJ J LOi,...,0J4 



(6) 



with x^2 for field-independent wave-function renormalization and 



u M 2(wi,w 2 ) = ^y^uiu) 2 + |*>(w 2 + u 2 ) + ^ = vh,wi) =v li 2(-u} 1 ,-w 2 ), (7) 
where the first two terms stand for the field-dependent, quadratic wavefunction renormalization, 

1- ,2-2 . l-o „2/-, 



2^ a O +2 Y n 2 l(M))> ( 8 ) 
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and the third term for the quartic self-interaction. Integrating by parts one finds 

- J^a,tq a ,tqb,tqb,t = J[q a ,tQa,tQb,tQb,t + 2Qa,tQa,tQb,tQb,t\, (9) 

so that one can perform the replacement io\ <^ + 2wiw 3 ) in the kernel of the last integral term in the right-hand 

side of Eq. (6) and write 

- l(lv^i 2 k 2 + l^g 2 (M)) = - l(ly^i 2 k 2 + l Y Am) 2 ) ( 10 ) 

= - Yp. 2 and Y n 2 = -2?>- 
Let us separate off the diagonal piece G -1 and the off-diagonal piece A of the matrix 

TW\$+w 2 = G- 1 +A, (11) 

where G (cA ,) i(di( y) = G<A,d<L+w',o with G u = [x^oj 2 - n^+ifi 2 }- 1 and 

S 2 f 

oqc,—u)Oqd,-u)' Juii,...,u4 



-w ^a;i,...,W4 

with 



= -2 / £_ w _ w / +Wl+W2j o[t^2(wi, W2) + V^(tjJ, w')]9a,wi3a,w 2 = -Bw'.w, 
u;i ,CJ2 

^(c,w),(d,a;') = _4 / tf-w-w'+wi+u^ofaj* 2 (^1, -w) + 1^2 (^2, -^'^cui 5d,w 2 = <5(d,«'),(c,a;)- (13) 

As to the next we Neumann-expand the inverse matrix in the right-hand side of the evolution equation (5), keeping 
the terms up to the second order in the off-diagonal piece, 

-hriG-i+A}- 1 = -^tr[G - GAG + GAG AG - + ...] 



Nh 



I G u + -( G£-A(o,w),(a ,-«) - 2 / G , ^A( CiW ) i ( d ^/)G aJ 'A( d: _ w /) ! ( C! _ w ) - + . . . (14) 

Jo; ^ J u ^ J lo,lo' 



It is worthwhile noticing here that for N = 1 the terms for field-dependent wavefunction renormalization given in 
Eqs. (8) or (10) merge into a single one, 

~ / S ull+ ... +UJ4 fl-T ^ui^q^q^q^ (15) 

with T^2 = y„2 — 3Y^2 . It is illuminating to give the graphical representation of the first few terms of the Neumann- 
expansion in the functional CSi equation (5) for the case N = 1, shown in Fig. 1. 

The details of the evaluation of the integrals in the right-hand side of Eq. (5) are given in Appendix A 1. In addition 
to the truncation of the Neumann-expansion we have expanded each term of the Neumann-series in powers of the 
Fourier-transform q a ^ of the field variable as well as in powers of its frequency index u and neglected the higher-order 
terms absent in Ansatz (6). In order to perform the expansion in the frequency (i.e. in the time-derivative of the field 
variable) we used the expansion of the propagator, 

G w+Q - G LJ -2x^u J Gla-[l-Ax^oj 2 G LJ }x^Gla 2 + 0(a 3 ). (16) 

Furthermore, we have to make use of identity (9) in order to transform all derivative terms into one of the forms present 
in Ansatz (6). The truncation the Neumann-expansion and that of the gradient expansion may lead to erroneous 
results in the IR scaling regime for sufficiently strong coupling g B and/or large number N of dimensions. The steps 
described above enable one to express the trace in the right-hand side of the evolution equation (5) in terms of the 
loop integrals 

I n ^ 2 ,nl 2 ) = hj A ^GZco s (17) 
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FIG. 1: Feynman diagrams for the vertices and the right-hand side of the evolution equation in CSi RG for the N = 1 
dimensional oscillator. 

and find the evolution equation 

-J 5w 1 +...+u i fi(^p^y IJ 2UJiL02 + 2^ 2 ^ 2W l + ^>^2g^jq a ^ 1 q a ^ 2 q b ^ 7i q b ^ i 

TN 1 - f 

= — 2"/i,o + 2 i N y» 2 -( N + 2)v]/ 2 ,o y w?«„ 1&;i 



4([^-(^ + 2)V]/ a , a -^^ 



(18) 



where the purely quartic terms Ti (i = 1, 2, 3) are given by Eqs. (A9)-(A11). In the limit A — > oo the loop-integrals 
can be taken analytically, and all expressed in terms of a few ones, I n _o with n = 1,2,3,4, (see Appendix A 2). 
Comparing the coefficients of the corresponding terms of the gradient expansion in both sides of Eq. (18), one arrives 
at the RG evolution equations for the couplings in the our approximation, called previously AWF, 



fl N T 



(19) 
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d^n 2 = -[Ny - (N + 2)F]J 2 , 2 + (iV + % J 2 , , (20) 

a M 2 ff = 12[iVy 2 + 2(JV + 2)yV + (TV + 8)V 2 ]/ 3 , 4 - 4 ff [(7V + 2)y - (N + 8)Y]/ 3 ,2 + ( N J 30 , (21) 

S^z = [JVy-(JV + 2)y]J 2 , , (22) 

0^ = 8x 2 [iVy 2 - 2(JV + 2)yY + (N + 2)V 2 ]J 5 , 6 - 8(jV + 2) V (y - y)J 5 , 4 + 2(iV ± 2)g V / 5 , 2 

-10x[Nf - 2(N + 2)yY + (N + 2)V 2 ]/ 4 , 4 + 2(JV + 2) 5 a;(y - Y)h. 2 - ^±^- x I ifi 

lo 

+4[-2y 2 — (N + 2)yY + (N + 2)Y 2 ]/ 3 , 2 + {N ^ (y + Y)I 3fi , (23) 

d^Y = 8x 2 [Nf - 2{N + 2)yY + (N + 6)Y 2 ]I 5 . 6 - ^-x 2 [(N + 2)y - {N + 6)Y]J 5 , 4 + ^±^x 2 / 5 , 2 

-\0x[Nf - 2(N + 2)yY + (N + 6)V 2 ]/ M + 2gx[{N + 2)y — (N + 6)?]7 4;2 - i^±3^ xh o 

lo 

+ [2iVy 2 - 8(JV + 3)yY + (6N + 40)Y 2 ]/ 3 , 2 - ^-[2y - (N + 7)Y]I 3fi . (24) 

For the sake of simplicity, we have suppressed the lower index /i 2 of the running couplings. The loop-integrals 
occurring in these equations are UV finite except of 7i ; o, so that the IR limits of all couplings can be determined 
for A — > oo with exception of the constant term of the 1PI effective action determining the energy of the ground 
state. The x-dependence of the right-hand sides of the evolution equations is partly explicit, partly implicit occurring 
via the propagators in the loop-integrals. The evolution of the field-independent wavefunction renormalization x is 
generated via the evolution of the field-dependent wavefunction renormalization parametrized by y and Y (c.f. Eq. 
(22)) because those start to evolve already at the UV scales due to the terms containing the loop-integrals I$ i2 and 
7 4> o in the right-hand sides of Eqs. (23) and (24) even for vanishing bare values ys = Yg = 0. The evolution equations 
for the LPA can easily be obtained by setting x = 1 and y — Y = identically It is worthwhile noticing that the 
analytic continuation of the RG equations to Euclidean space via the replacement —ifi 2 — > A for real A decreasing 
from Xb ~ A 2 towards the IR limit A — > is rather suitable for the numerical treatment since the Euclidean RG flow 
does not develop imaginary parts for the couplings. 

For N = 1 only the single coupling T = y — 3Y for field-dependent wavefunction renormalization occurs and the 
RG equations (19)-(24) reduce to the following ones, 

9^1 = \h,o, (25) 

d^n 2 = -ri 2 , 2 + \ghfi, (26) 

d^g = 12T 2 / 3i4 -12 5 T/ 3i2 + 3 5 2 / 3 ,o, (27) 
d^x = TJ 2 , , (28) 
d^T = -16.T 2 T 2 / 5 , 6 + l&gx 2 Th A - 4 5 V/ 5 ,2 + 20T 2 a;/ M - 12 ff :rT/4,2 + g 2 xI 4 ,o - 14T 2 / 3j2 + 5.gT/ 3 , . (29) 



We see now for N — 1 that the parameter x starts to evolve in the UV scaling regime because of the parameter T 
does due to the nonvanishing terms — ^g 2 x 2 I^^ 2 + g xlijo in the right-hand side of Eq. (29). 
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B. Solution of the CSi RG equations 

1. IMA 



The IMA corresponds to the replacement of the running couplings by their bare values D, B — Uq, gs, xb = 1, 
%Ib = Yb = into the right-hand sides of the evolution equations (25)-(29). Making use of the loop integrals 
of Appendix A 2 in the limit A — > oo and the analytic continuation to the Euclidean space via the replacement 
A = — i/x 2 , one can easily integrate Eqs. (25)-(29) from the initial scale A B — > oo down to the gliding scale A. The 
UV divergence of the ground-state energy can be removed by the choice -fs = \NH\uj\ 1 + As] 2 . The UV scaling laws 
obtained in that manner can then be extrapolated to the IR limit A — > 0: 



1 + 



4(N + 2)£ 



go = gB 



2{N + 8)£ 



1 4(A + 2)u o e - 
x a = 1, hy = , HY a 



4{N + 6)w C 2 



(30) 



where the dimensionless parameter £ = of the Rayleigh-Schrodinger perturbation expansion [39] has been in- 
troduced. We see that the one-loop corrections of the couplings depend linearly on the dimension N. The en- 
ergy of the ground state Eq = 70 = ^NhujQ is just the sum of the zero-point energies of the number N of in- 
dependent harmonic oscillators in this approximation, the one- loop result. The energy gap between the first ex- 
cited state an d the g round state, i.e. in held-theoric terms the minimum energy of the one-particle excitations, 
Ei — Eq = hy/Q^/xQ = + 4(N + 2)(£/3) increases with increasing bare coupling strength gs and number of 

dimension N. The extrapolated IR value of the quartic coupling go decreases with its increasing gs and N in the 
IMA. Although there occurs a held-dependent wavefunction renormalization of the order O (g%) in the IMA, the 
field-independent wavefunction renormalization is absent since it is a two-loop effect in the lowest nonvanishing order 
of the loop expansion. The beta-functions do not change sign, so that one can conclude that the extrapolation of 
the UV scaling laws to the IR region seems to predict only a single phase of the quantum oscillator independently of 
its dimension N. However, restricting ourselves to that approximation, we would get in trouble with the change of 
sign of 50 for sufficiently large gB and N, as well as with a continuation of these formulae to negative values of ujq 
corresponding to a symmetric double-well bare potential. 



2. Numerical integration of the CSi RG equations 

The numerical integration of the RG equations (20)- (24) obtained in the CSi RG scheme in the AWF and analytically 
continued to the Euclidean space has been performed by means of a fourth-order Runge-Kutta algorithm. In Figs. 
2 and 3 the RG flow obtained in AWF is compared to RG flows obtained in the IMA and the LPA. (The latter is 
obtained by setting x = 1, y — Y = identically.) One can see, that the flows in the various approximations are 
qualitatively the same, i.e. after the UV scaling region both the dimensionful frequency fi 2 and the quartic coupling 
g saturate at constant values in the IR region. The UV scaling is described properly by all approximations, but they 
give different values in the IR limit. The IMA corresponding to the one- loop approximation gives the largest value for 
fi 2 and the smallest one for g. More reliable results are obtained in the LPA and AWF which resum infinitely many 
Feynman diagrams. Both of the Figs. 2 and 3 clearly show that the inclusion of the wavefunction renormalization 
modifies only slightly the evolution as compared to the LPA. This is a consequence of the fact that the pole of the 
propagator, corresponding to the first excited state of the oscillator dominates in its Lehman-expansion, as discussed 
in [27]. 

The flows in the AWF have common features that can be recognized from Figs. 2-5. Since the couplings y and Y 
scale quite similarly, we have plotted their combination T into which they merge in the case of the one-dimensional 
oscillator. After an UV scaling region extending roughly over three orders of magnitude of the scale parameter A, the 
flow saturates in the IR region for any of the dimensionful couplings, so that all couplings are IR relevant. Namely, 
they tend to nonvanishing constant values in the IR limit and only a single phase is detected. The latter corresponds 
to the general expectation because there is no room for spontaneous symmetry breaking in the 0+1 dimensional 
quantum system. 

It has been found that the flows for any of the couplings do not alter qualitatively with increasing N. That is 
illustrated in Figs. 4 and 5 for the parameters of the wavefunction renormalization. The A^-dependence of the various 
couplings in the IR limit is depicted in Fig. 6 for the various parameters. In this respect the IMA turned out to be 
misleading. According to the IMA the quantum fluctuations in the IR limit result in additive contributions to the 
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FIG. 2: The flow of the frequency parameter Q. 2 in different approximations, for Q,% = 0.00102, qb = 0.0001, xb = 1, 
VB = Yb = 0, and JV = 1. 
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FIG. 3: The flow of the quartic coupling g in different approximations, for 0,% = 0.00102, qb = 0.00001, xb = 1, J/s = Yb — 0, 
and JV = 1. 



bare values of the couplings, which are linear in the dimension N, but the AWF gives completely different behavior. 
We see in Fig. 6 that the parameter il§ shows up a powcrlike dependence on N for sufficiently large dimensions, 
fig oc N bn with bn w 0.58. In field-theoric terms the oscillator gets more 'massive' with increasing dimension N. 
The coupling g depends rather weakly on N having a maximum at around N w 30. The slow decrease for large N 
values may be the consequence of the truncations used. The neglection of higher-order monomials in the gradient 
expansion can cause a failure when one or both of the dimension TV and the bare quartic coupling g B are too large. 
It is easy to notice from the explicit form of the RG evolution equations in the IMA that the loop-expansion goes 
effectively with Ngsh for large N, so that the higher-order terms of the local potential can grow up to the same order 
of magnitude as the quadratic and quartic terms are when the flow reaches the border of the UV scaling region and 
that may result in their saturation in the IR scaling region at values of the same order of magnitude, too. In Refs. 
[26, 27] the authors take with the higher-order monoms of the local potential, as well. They have found that the 
renormalized coupling go increases strictly monotonically with increasing bare parameter gs for the one-dimensional 
oscillator. Similarly, g should increase with the dimension N of the oscillator also strictly monotonically, because 
the the loop-expansion goes effectively with Ng B for large N. As one can see in the inset of Fig. 6, the deviation 
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FIG. 4: The scale dependence of x for Q% — 0.001, qb = 0.01, and y = Y = 0. The order of the curves from the top correspond 
to the various values N = 2, 10, 20, 50, 100, 200, respectively. 
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FIG. 5: The scale dependence of T for 0,% = 0.001, qb = 0.01, and y — Y — 0. The order of the curves from the top correspond 
to the various values N = 1, 10, 20, 50, 100, 200, respectively. 



xo — 1 of the field-independent wavefunction renormalization parameter xq from its bare value as well as the coupling 
To of the field-dependent wavefunction renormalization go to zero according to the power laws, xq — 1 ~ N~ bx and 
To ~ N~ hr with b x w 0.55 and by w 0.45, respectively for increasing N. Therefore the wavefunction renormalization 
effects seem to be washed out in the limit N — > oo because of the increment of the mass gap. 

The dependence of the IR parameters on the bare ones can be characterized by the sensitivity matrix, which is 
defined by 



81,521 ' ' d S2 ' 



(31) 



where s\ — VL 2 , g, x, T are the running couplings and s 2 = £1%, gs, %b, Tg, and N are the bare parameters 
and the dimension, respectively. The elements of the sensitivity matrix integrate the infinitesimal changes of the 
infinitesimal blocking steps during the RG flow when the ratio r = A/As of the moving scale to the UV scale glides 
from r = 1 to the IR limit r = 0. Therefore, their values in the IR limit, S SliS2 (0) as characteristics of the global 
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N 

FIG. 6: The iV-dependence of the IR values fig (full line) and g (dashed line) for Q.% = 0.001, g B = 0.01, and y = Y = 0. The 
IR values of the wavefunction renormalization parameters xo (full line) and To (dashed line) are shown vs. N in the inset. 
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FIG. 7: The flow of the sensitivity matrix element Sq2 N for fl B — 0.001, qb = 0.00001, and ys = Yb = 0. The curves from 
above correspond to N = 1, 2, 5, 10, 20, 50, 100, 200, respectively. 



RG flow may show up singular dependences on the bare parameters, which may reveal themselves in qualitatively 
different scale-dependences S SltS2 (r) along the various RG trajectories started in various regions of the space of the 
bare couplings [41]. Since the model in our case has a single symmetric phase, no such singularities are expected to 
occur in the sensitivity matrix. This is demonstrated in Figs. 7 and 8, where the matrix elements S^.n an d S^ 2 ,^ 
are shown to vary monotonically with decreasing scale A. The flow of these matrix elements do not alter its character 
with increasing N (illustrated in Fig. 7). The UV and the IR scaling regions are again clearly distinguishable and the 
flow saturates in the IR region. The larger is the dimension N, the weaker is the sensitivity of the IR parameter £Iq 
to the same small increment A7V of the dimension. 

Quite similarly, the character of the scale-dependences of Sn2^ N (r) and Sq2 .n 2 B ( r ) do not alter starting the RG 
flow with various bare values Q 2 B (see Fig. 8), including also its negative values corresponding to a double- well bare 
potential. This is another signal that the anharmonic oscillator exhibits only a single, symmetric phase. As it was 
thoroughly discussed in [26, 27] the large-amplitude, low-momentum quantum fluctuations fill up the minima of the 
double- well potential and one ends up with a convex potential and a unique ground state in the IR limit. It was also 
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FIG. 8: The flow of the sensitivity matrix element S n 2 ^ for N = 50 and gs = 0.00001. The curves from below correspond 
to f2| = -0.001, - 0.0001, 0.001, O.Ol, respectively. 
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FIG. 9: The iV-dependence of the IR values f2§ vs. the bare parameter Q%. The curves from below correspond to N 
1, 10, 20, 50, 100, 200, respectively. 



discussed that for sufficiently small bare quartic coupling the RG evolution equations show up singularities and one 
has to take into account time-dependent saddle point configurations for evaluating the path integral, that has been 
achieved by turning to an instanton gas approximation. In Fig. 9 we illustrate that the RG flow started with negative 
parameter values Q B < ends up with some positive values in the IR limit. We note that in our case the evolution 
gets wrong for too large negative UV values of Q 2 B due to the strong truncation of the local potential in the ansatz (6). 
As N is increased we get larger and larger IR frequency values Qq according to the power-law established previously. 

Closing the discussion of the numerical results obtained in the CSi RG scheme, we can make an estimate of the 
spatial extension of the ground-state wavefunction. As shown in the Appendix A 4, the variance of the coordinate 
operator q can be estimated by making use of the Thomas- Reiche-Kuhn sum rule [40] and the dominance of the first 
pole in the Lehmann-expansion of the single-particle propagator [27]. With the help of the scaling of the renormalized 
coupling with increasing dimension N established above, we get from (A37) that the variance scales as 

( | £ 2 |0) « - ~ N 1 -^ ~ N - 71 . (32) 
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This means that with increasing dimension N the quantum oscillator becomes more and more delocalized because of 
the increasing the number of degrees of freedom faster than that of the 'mass', ie. the frequency fl of the oscillator. 



III. ISG RG FOR THE ONE-DIMENSIONAL OSCILLATOR 



A. Evolution equations 

Further on we shall restrict ourselves to the study of the one-dimensional anharmonic oscillator, since the scaling laws 
for the TV-dimensional anharmonic oscillator do not show up any significant qualitative modification in dependence 
on the number of dimensions N according to our findings in the framework of the CSi RG scheme. Now we choose an 
IS RG scheme in that the bare coupling gs of the anharmonic quartic interaction term acts as control parameter and 
it is evolved from zero towards a final positive value. For the sake of simplicity, below we shall suppress its index and 
denote the control parameter of the evolution simply by g. Here we apply the ISg RG scheme to the so-called second 
Legendre-transform, the Legendre-transform of the generating functional W[J] of the connected Green- functions of 
the bilocal operator qtqt> defined via the path integral 

< iww> - j^mw-m-^-^v^ pa, 

with the symmetric external bi- local source J t t > = J t > t [35]. We have introduced the infinitesimal imaginary piece of 
the quadratic term of the potential by the replacement Wq — > uj 2 — ie in order to make the path integral convergent. 
In the presence of the external source J the propagator G is given as 

W$ = ™ = (T( qtqt ,)) = G M „ = = G _,. (34 ) 

The vanishing source J = corresponds to the ground state with the propagator (T(q t qt>)) j=o — G gr t,v ■ The second 
Legendre-transform, the 2PI effective action T[G] is defined as 

r[G]--^[J] + itr(G tr - J), (35) 

and the relations 

ST[G] 1 8T[G] 1 

6G t , t , ~ 2^*'' SG^ 2 J W 

hold. The symmetry of the source J implies the symmetry Gt,t> — Gt>.t of the propagator. The second functional 
derivatives satisfy the identity 

I = W {2) :T {2 \ = r (2)-i/ = r (2)-i :J ( 37 ) 

in shorthand notations with A : B = f , A(^( ttl ^B( tt ,)^\ l = J u , ^(..^.(^.^^^(-aj.-w'). (...)• For later convenience 
we introduce the matrices /, L, and 1, 

I(t 3 ,U),(ti,t 2 ) ~ 2^* 1 '* 3 ^* 2 '* 4 + ^*i.*4^*2 ,t 3 )i ^(w 3 ,w4),(«i,ws) = 2(^1+^3^^2+^4,0 + ^1+^4,0^2+^3,0): 

^(t3,t 4 ),(tl,t 2 ) — ^t 1 ,t 2 8t 1 ,t 3 8t u t 4 : ^(U3,U 4 ),(U1,U2) — ^Wl+W 2 +W3+W4,0l 

l(t 3 ,*4),(tl,*2) = ^l,t 3 ^t 2 ,t4! 1(W 3 ,W 4 ),(W1,W2) = ^Wl+W3,0^W 2 +W4,0- (38) 

The matrix I acts over the two-particle (in field-theoric sense) subspace as the projector on the symmetrical subspace 
(that of the two-particle states being symmetric under the exchange of the particles) and plays the role of the identity 
operator in that subspace and L : I — L holds, while 1 represents the usual identity operator over the entire two- 
particle subspace. The inverse A -1 of an arbitrary two-particle matrix A satisfies A^ 1 : A = 1, while for the inverse 
A^ 11 = A^ 1 : I of the same matrix in the symmetrical subspace A~ 1! : A = I holds. Further on let us denote the 
trace of an arbitrary two-particle matrix by TrA = j t t , ^4( t ,t'),(t,t') = f UiU , ^(«,«'), (-«,-«')• 

The bare coupling g is turned on gradually from zero to some finite value in order to control the evolution of the 
2PI effective action satisfying the following RG evolution equation, 



d g T[G] = -d g W[J\ = -^e-* w W J Vq ^qfe^^<i 2 -^a-^9 2 -^)+^i-'-9 



24 W 



l -Tr(L: W™) + ( ^ ) W^:L:W^ 



2 



i_h 

247 



Tr(L : r^" 11 ) + l -G : L : G 



(39) 
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Separating off the trivial term controlling the RG evolution, one can introduce the reduced 2PI effective action T[G] 
via the relation 



T[G] = f [G] + ±G:L:G 



and recast Eq. (39) in the form 



(40) 



(41) 



B. Solution of the evolution equation 

1. IMA 



In order to make a guess on the functional form of the reduced 2PI effective action f [G] for the linear quartic 
oscillator, we evaluate it first in the tree-approximation, and solve the RG equation in the IMA. The tree-level 
approximation corresponds to the vanishing value of the control parameter, g = 0, ie. the case of the linear harmonic 
oscillator, for which the Gaussian path integral can be evaluated in a straightforward manner, yielding 



■it) 

W trcc [J] = ^trln[(w 2 - u $+ie)I + J] 



(42) 



and 



r tro °[G] = —tr\nG--tY[{uj 2 -ul + ie)I-G]+ const. = f troo [G] 



(43) 



up to an UV divergent additive constant (being independent of the bare parameters). The necessary condition of 
the extremum of the 2PI effective action in the tree-approximation, ST IMA [G]/SG = provides the propagator in the 
ground state of the linear harmonic oscillator, 



as expected. The inverse of the noninteracting two-particle propagator is given via 



/-jtree— II 

y (— U>3,— U>4),(— Ul,U2) 



[G] 



=,tree(2) 

(— Ul ,- W2),(— U13,— UI4.) 



[G] 



ih 



s-1— 1 f^ -1 1 



(44) 



(45) 



Hence one finds 



(46) 



that we write as g trcc = l^QG) 1 introducing the shorthand notations {AB)^ c d ^ — ^(A c . a Bd,b + A^ a B Ci b), and 
for later use a :: (5 — J a d (X[...] l [(c,d),(a,b)]P[(-c,-d),(-a,-b)],[...] where Latin indices stand for frequency indices. For 
later use we also evaluate the first and second functional derivatives of g tlee [G] , 



sg 



[G] il 



SG- 



+ S Ui+u>0 d U2+ult0 G U3tUl + (w 4 UJ 3 )} + (LJ <^> Lj') , (47) 



s 2 gy cc w AG] 
5G SG 



U.—W u*->—u" —u'" 



Ih 



[^W3+W 1 0^W1+W',0^W4+W",0^W 2 +W'",0 + ^W4+W,0^W2+W',0^W3+W",0^W1+W"',0 



+ (uj 4 w 3 )] + (w u') ] + (uj" ^ uo 



(48) 



In order to obtain the 2PI effective action in the IMA, we have to insert the tree- level expression r tree [G] into the 
right-hand side of the evolution equation (41) and integrate it over the control parameter from the initial value g = 
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to some finite value g. Then we get 



1 h 



r 1MA [G] = T trcc [G] + — j I dgT: 



r tree(2)j G ] + ^ L | . £ 



-II 



I :ln(l + ^L:£ tree [G] 



= f tree [G] + f gU^fX) iv(i:^ r;] 



2 ^ n V 12 

n=l v 



(49) 



The inverse propagator of the ground state in the IMA is given through the necessary condition of the extremum of 
the 2PI effective action r IMA [G], ie. <5r IMA [G]/<5G = 0, that yields the equation for the propagator of the ground 
state in the IMA, G™ A , 



G 



IMA -1 

gr 



Gg r °° _1 + ^Tr 



G=G™ A 



i "y u . piMA\ 
+ 6h {L - G * c 1 



(50) 



ie. 



/-yIMA —1 /ttree— 1 

gr uj,uj' gr 



ig_ 
2h. 



«/wi,W2 

~l~ -r Qfc2 / ^+Wi+W3+CJ4,0^'+W2+CJ5+^6,oG t j 1 . t j 2 G W 3. W5 G W4!W6 + O (tj ) 



G=G™ A 



(51) 



when the terms up to the order O (g 2 ) are only kept. Multiplying both sides of this equation by G™ A from the 
left and G* r ee from the right, one would obtain the Schwinger- Dyson equation in the IMA for the propagator of the 
ground state, but for our later purposes it is more convenient to retain the present form of the equation. The solution 
of Eq. (51), the propagator of the ground state can be used to determine the energy of the ground state, 



£o= T lim ^r IMA [G^ A ] 



(52) 



and the 2PI 4-point vertex function, 



[2] IMA _ £IMA-1I 

'(u>3,U>4),(wi,W 2 ) y (u 3 ,U 4 ),(ui,U2) 



■IMAi 



-'tree— II 



-yIMAi 



St 

— ^°Wl+O)2+W3+W4,0 



(w 3l U4),(ui,u 2 )l °gr J 

~ / [^wi+ws+wJ+w^O^+Wi+Wa+^O + ( Wl ^ U>2 ^ G e^ t[,u' 2 G ™ t' 3 ,w 4 + (.9 3 ) 

JuJi,a;2,a;3,aj 4 



(53) 



with g IMA - u [G] =r IMA P)[G]. 

The evolution equation (41) allows for a simple diagrammatic picture, shown in Fig. 10. The single vertex without 
the factor g (empty circle) present in each of the diagrams for d g f occurs due to the partial derivative with respect to g. 
The result of the IMA is the rcsummation of the ring diagrams shown for T[G] — r trco [G] in the fourth and fifth lines in 
Fig. 10. In the fourth line we depicted the expansion of the 2PI effective action in powers of the tree-level two-particle 
propagator Q trec , each diagram consists of double-line sections for Q tree, s joined via the vertex functions 7 PI trcc to a 
closed loop. As one can see from the expansion in powers of the propagator G shown in the fifth line, the 2PI vertex 
function is taken at the tree-level, 7^ trcc = gL/ 12. The proper self-energy insertion £ IMA — G _1 IMA — G _1 tree and 
the 2PI 4-point vertex function 7 PI IMA in the IMA are depicted in Fig. 11. 



2. Expansion in the 2-particle channel (E2PC) 

The functional dependence of the reduced 2PI effective action f IMA [G] on the tree- level 2-particle propagator g tlee 
in the IMA brings one to the idea to look for a better approximation of the effective action as an expansion in powers 
of the tree- level 2-particle propagator g tTCC = 2^(GG) J making the Ansatz 

f E2PC [G] = r troo [G] + ^j [2] :: (GG) 1 + ^(GG) 1 :: 7 [4] :: (GG) 1 + ..., (54) 



16 



e= ~ +X -n 




FIG. 10: Feynman diagrams for the RG evolution of the 2PI effective action in the IMA. 




FIG. 11: Feynman diagrams for the proper self-energy insertion S = G 1 IMA — G 1 trcc and the 2PI vertex function 7 [2 ' IMA . 
where the symmetry properties 

[2] [2] _ [2] 

!{a,b),(c,d) ~ i(b,a),(c,d) ~ '(c,d),(o,6) ' 

[4] _ [4] _ [4] _ [4] 

'[(o,6),(c,d)],[(o',6').(c',<i')] ~~ ~[(M),(c,d)],[(a',6'),(c',<i')] ~~ ~[(c,d),(a,6)],[(a',6'),(c',d')] — '[(a' ,b'),(c' ,d')U(<*,b),(c,d)] y 00 ) 

are required (Latin indices stand for frequency indices) for the 2PI vertex functions, 7^ and 7M. We shall solve 
the evolution equation (41) inserting the Ansatz (54) into it, expanding its right-hand side in Neumann-series, and 
keeping the terms up to the quartic ones in G, ie. those of the order O ([(GG) 7 ] 2 ) in its both sides. The approximation 
with the E2PC goes beyond the IMA when the inverse of the 2-particle propagator, r( 2 ) has been replaced by 
r trcc(2) = gtree-i[ G ] in thc right-hand side of the RG equation (41). Now the Ansatz (54) in powers of g trcc [G] 
implies 

f E2PC (2) [G] = Qtree-1 [G] + -[2] + q (( GG )/) . ( 56 ) 

Having expanded the inverse matrix, ic. 2-particle propagator in the right-hand side of (41) in Neumann-series as 

g E2PC [G] = r E2P C(2) -1[ G ] = ^fE2PC (2) [G] + 1 = gtr CC[G] _ g trcc [G] . ^2] . g trcc [G] + Q ([( G G) J ] 3 ) (57) 

with the 2PI 4-point vertex function, 

7 [2] = p(2) [G] _ gtree-1 [G] = - [2] + Jr ( 58 ) 
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we obtained the following coupled set of RG evolution equations for the 2PI 4-point and 8-point vertex functions, 



Q =J2] = 1 J 

U 9 '(W1,W2),(W3,W 4 ) g-^(w3,W4),( w l> w 2)' 



aV 4] = -- 

9 '[(wi,u 2 ),(w3,W4)],[(a;;,w^),(w^w^)] ^ 



+-k(w3,w 4 ),K,^)U : T [2 ' 1 ^)(wi,W2),(wl,w 2 ) ) + (W1,W 2 ) ^ (w 3 ,W4 



(59) 



Those can be integrated straightforwardly for the initial conditions J^l — 0, 7^0 = 0, and one gets 



-[2] ±J 
T(W1,W 2 ),(W3,W4) ~~ g lv (wi,W2),(w3,W4)' 

„2 



'T[((l;i,W2),(w3,W4)],[(w'i,w 2 ),(^,W4)] — 4^ [^wi +w 2 +o4+w 2 ,0^W3+W4+w 3 + w 4 >° + ^3 +"4 +W3+W4 ,0<Li +w 2 +u' 3 +w 4 ,o] (60) 

and hence 7^ = |L. Since we have kept only the contribution of O {g) to the reduced 2PI 4-point vertex function 7 PI, 
the expansion (54) of the functional f E2PC [G] in powers of the free 2-particle propagator Q tree <~ (GG) 1 corresponds 
to its expansion in powers of g. The 2PI effective action takes now the form 

r [G] r [G] -\- '— I ^! + ,0^Wi ,W2 ,W4 

° J U>\ ,...,UJ4 

— Zsfi / ^i+w2+w;+w 2 > ^3+w4+wj+w;,oG llI1 , U3 G l j 2 , [ j 4 G IJ ' iill ,jG lJ ^^ + (G 6 ) . (61) 

In this approximation based on the E2PC there occurs an additional contribution of the order O (<? 2 ) in the 2PI 
effective action which was not included in the IMA: in the right-hand side of Eq. (61) there occurs an additional 
factor 3 in the term of the order O (g 2 ) as compared to the corresponding term in the IMA. This is because of the 
improvement of the 2-particle propagator, 

g E2PC [G] = g trcc [G] + |g trco [G] : L : G trcc [G} + O (g 2 ) (62) 

(c.f. Eq. (57)) as compared to the free one used in the IMA. Putting it in another way, in the approximation base on 
the E2PC we use the 2PI 4-point vertex function 7^ = 7 PI + 7 PI tro ° = [g/A)L instead of 7 PI troc = (g/12)L used in 
the IMA. The diagrammatic form of RG the evolution equation remains the same as that in the IMA (see Fig. 10) 
except of replacing the vertices 7^ trcc by 7^] of the present approximation. 

The necessary condition of the extremum of the 2PI effective action, <5r E2PC [G]/<5G = provides again the 
Schwinger-Dyson equation for the propagator of the ground state, 



pE2PC -1 _ ^trcc-l 



IE. 

2ft 



9 f 

I ^wi+o;3+cj^+a; 2 . 0^2+^4+^3+^4, 0^^3,^4^Wi,W3^ ,f ^ 



/ U}{,..., W 4 ,W 3 ,W4 

and a straightforward but lengthy calculation yields 

[2] E2PC _ r E2PC -II r/-<E2PCi _ r troo -II r r .E2PC 



(63) 

G=G E2PC 



'(W3,W4),(W1,W2) (W3,W4),(W1,W2) ^ S r ' (W3 1^4) , (lUl ,U>2) ^ S r ' 

= 2^1+"2+W3+W4,0 - -jTT / [^Wl+W3+Wi+ W 3.°^2+W4+W 2 + W 4. + ( Wl ^ W 2 )] [G w J ^ G w ^ J G = G E2PC (64) 



for the value of the 2PI 4-point vertex in the ground state. The O (g 2 ) term of the propagator as well as that of the 
2PI vertex function exhibit the additional factor 3 as compared to the corresponding terms in the IMA. 



C. Determination of the renormalized couplings 



We shall obtain the renormalized quantities in the IS RG scheme in the IMA and in the truncated E2PC, as 
described previously in Sections IIIB1 and IIIB2, respectively, in both cases keeping the terms up to the second 
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order O (g 2 ) of the bare coupling g that has been used as the control parameter. We shall make use of the diagonal 
form of the propagator, G UtU! i = G w £ WjW '. First, we determine the propagator G gr u of the ground state by solving 
the Schwinger-Dyson equation with the precision O (g 2 ) , and afterwards the energy of the ground state via 



£0= lim -r[G gr ], 

1 —too 1 



(65) 



and the renormalized coupling go from the 2PI 4-point vertex function taken at the symmetric point u>i = (i = 
1,2,3,4) via 



lim - [2] 



(66) 



Finally, we compare these results with those obtained at the same order of the perturbation expansion in the CS 
RG scheme given in Appendix A3, in the WH RG scheme given in Appendix B2b, as well as by the results of the 
well-known Rayleigh-Schrodinger perturbation expansion in quantum mechanics. 

First, let us present the results obtained in the framework of the IS RG in the IMA. For diagonal propagator the 
Schwinger-Dyson equation (51) can be rewritten as 



G 



IMA -1 

gr bJ 



Gtree —1 , 
gr u + 



2ft 



/ G Ul + 2 / S LJ+UJl+L j 2+LU3 fiG UJl G UJ2 G UJ3 + O (g 3 ) 



G = G™ A 



(67) 



Let us insert into Eq. (67) the perturbation expansion 

G- 1 = G 1 +gA + g 2 B + O(g 3 ), G = G - gG ■ A - G - g 2 G ■ C ■ G + O (g 3 ) 



(68) 



where A and B are symmetric matrices to be determined, C = B — A ■ Go ■ A and we wrote Go = G g ^. cc for the 
sake of simplicity. Comparing the terms of identical powers of the coupling g in both sides of Eq. (67), one finds 
A u ,uj' — A^Scj+uj' ,o an d B^.u' = BuSu+ur .o with 



— ^ J Go Wl , B u — ^2 J Go W1 J Go ui 2 + Yg^2 J ^° "i+w 2 -wGo Ul Go 



The loop integrals can be taken explicitly making use of the residuum theorem, 

[ Go lu = - — , f Go w = — g , / Go m+uj-uGo wjGo ui 2 
Ju, Zu > 4cj J Ul , U2 



m 3 



Au: 2 w 2 - 9w 2 + 18ze 



(69) 



(70) 



Expanding the last loop integral in powers of uo 2 — ojq and keeping the terms up to the order O (oj 2 — uiq), one finds 
for the inverse of the propagator of the ground state 



G 



■IMA -1 



= --(x lu 2 - SIq + ie) 



with 



O 2 - , ? 
"o — w o 



0(g 3 ), x = l + -e+O(g 3 ) 



(71) 



(72) 



in terms of the dimensionless parameter £ = yf^a , introduced previously. For diagonal propagator G the 2PI effective 
action given by Eqs. (40) and (49) and the 2PI 4-point vertex function (53) reduce to 



r 1MA [G] = r tree [G] + ^T 



G u 



144ft/ 



-T 



I 



G U i 1 G U2 G U j 3 G U j 1 -\- U ;2+i03 ~t~ ® {G 



and 



[2] IMA 

^(W3,W4),(W1,W2) 



g w 

4 24ft 



/ G UJ (G- LJl - 

J U) 



G_ 



Ul2—U> 3 —U>, 



>U>1+U!2 +U>3 +U>A ,0 J 



(73) 



(74) 
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respectively. Inserting G — G™ A into the expression (73) one finds for the energy of the ground state 



E = -hu>o 



l+e-fe 2 +0(.9 3 ) 



According to Eqs. (66) and (74), one obtains for the renormalized quartic coupling, 



.90 = 9- i ^j^Gl„ + 0{ g *)=g 



l-^ + 0(.9 2 ) 



(75) 



(76) 



At second, we turn to the approximation based on the E2PC. Making use of Eqs. (61), (63), and (64), one can 
proceed like in the case of the IMA and obtain results with slightly different factors for the second-order terms, 



O 2 - , ? 



71 



(77) 



E = -huj 



.9o = .g(i-4£ + 0(£ 2 )). 



(78) 



(79) 



Finally, in Table I we compare the results obtained in the various approximation schemes. The various observables 
are given keeping the terms up to the order O (g%) of the bare coupling. We compare here the renormalized values, 
ie. the IR limits of them obtained in the CSi RG scheme in second order of the perturbation expansion as given in 
Appendix A3, the renormalized values obtained in the ISg RG scheme in IMA and in the approximation based on 
the E2PC when the evolution has been ended up and the control parameter reached the same bare value gs that used 
in the CSi RG scheme. We also present the results obtained in the framework of the WH RG scheme as outlined in 
Appendix B2b, and the results of the Raylcigh-Schrodinger perturbation expansion (RSPE) [39]. In the RSPE one 
finds for the energy of the n-th energy niveau (n = 0, 1, 2, . . .) of the one-dimensional anharmonic oscillator with the 
bare Hamiltonian H — ^p 2 + \u 2 x 2 + 



24 ' 



E,, 



n - 



,(2) 



(80) 



with 



,(2) 



(1 + 2n + 2n 2 )£, 

~7^ 2 Y, -T^—W 2n - - l)<W-2 + n(n - l)(n - 2)(n - 3)<5„,,„_ 4 



18 ~ " w - n~ 

n 7=n 

+4(2n + 3) 2 (n + l)(n + 2)5 n ^ n+2 + (n + l)(n + 2)(n + 3)(n + 4)5 n / )n+4 ]. 
This yields for the first two energy levels 



(81) 



En — 



Wi+£-^ 2 +o(£ 3 ; 



Ei - ^o(3 + 5C-|e 2 + 0(e 3 )) 



(82) 



and 



(hn ) 2 - {Ex - E ) 2 = (hujofil + 4£ - 12^ 2 ). 



(83) 

The same expression for the energy E of the ground-state has been found in Rcf. [37] by path-integral method. 
There are no asymptotically free states of the particle moving in the quartic anharmonic oscillator potential, therefore 
the definition of scattering amplitudes and through those a renormalized quartic coupling is problematic in the RSPE 
approach. 

All the RG schemes reproduce the first-order corrections in accord with the RSPE. This is because these corrections 
are one-loop corrections according to the loop expansion and the essentially one-loop RG equations describe the 
quantum effects at the one-loop order exactly, independently of the approximation schemes additionally used. We 
see that all the RG schemes have it common that the second-order corrections deepen the energy levels, tend to 
decrease the quartic coupling and result in a field- independent wavefunction renormalization x > 1. Nevertheless, 
the truncation of the perturbation expansion causes a strong dependence of the quantitative results on the RG scheme 
used. 
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Observable 


ISg RG IMA 


ISg RG E2PC 


CSi RG 


WH RG 


RSPE 


e q - ±hojo - (, 


55 £-2 
9 


16 t-2 
3 ? 


71 t2 
30? 


-18C 2 






55 t'2. 
6 * 


71 f2 
8 S 




-32£ 2 




S0-9B 
9B 








-6£ 





xo — 1 




if 2 


2 t2 
3? 









TABLE I: Comparison of the observables obtained in the framework of various RG schemes. The expansion in the dimensionless 
parameter £ — -^4r is truncated keeping the terms up to the second order. 



IV. CONCLUSIONS 

The role of wavefunction renormalization and dependence of the observables on the RG scheme in the second order 
perturbation expansion has been investigated on the example of the anharmonic oscillator. The functional CSi RG 
scheme with the imaginary mass and the ISg RG with the bare quartic coupling as control parameters have been 
applied to the 1PI and 2PI effective actions, respectively. The CSi RG evolution equation derived for the iV-dimcnsional 
0(N) symmetric oscillator has been solved by making a quartic Ansatz for the 1PI effective action including terms for 
field-dependent wavefunction renormalization, since that is the only way to generate field- independent wavefunction 
renormalization by means of a one-loop evolution equation. The evolution equation has been solved keeping the 
terms of the Neumann-expansion of the inverse matrix in the right-hand side of the evolution equation only up to the 
second order in the renormalized propagator and those up to the quartic ones in the oscillator coordinate. It has been 
shown that the O(N) symmetric anharmonic oscillator exhibits only a single phase, independently of its dimension 
N. Furthermore, the inclusion of next-to-leading order of the gradient expansion, the wavefunction renormalization 
does not modify qualitatively neither the IR scaling laws nor the phase structure. It has also been established that 
the effect of wavefunction renormalization decreases with increasing dimension N of the oscillator according to a 
power- law. 

The ISg RG has been applied to the 2PI effective action, in order to obtain a one-looplike RG evolution equation 
again. Our study was restricted to the one-dimensional oscillator in that case. An approximation scheme based on 
the E2PC, going beyond the IMA, has been put forward that corresponds the expansion of the 2PI effective action 
into the powers of the free 2-particle propagator g trec . The RG evolution equation for the 2PI effective action has 
been rewritten as a coupled set of evolution equations for the 2PI 4-point and 8-point vertex functions by keeping 
the terms up to those of the order O ((£ trco ) 4 ) in the effective action and those of the order O ((G tTCC ) 2 ) in the 
Neumann-expansion of the inverse of the full-particle propagator. The solutions for the vertex functions have been 
inserted into the Schwinger-Dyson equation for the propagator of the ground state that has been solved by making an 
Ansatz for the latter assuming the dominance of a single pole. The advantages of that ISg RG method as compared 
to the CSi RG scheme are that it enables one to read off the field-independent wavefunction renormalization directly 
and it provides an UV finite result for the ground-state energy. 

Finally we compared the results obtained for the ground-state energy, the 'mass' (in field-theoric term), the renor- 
malized quartic coupling, and the field-independent wavefunction renormalization in the various RG schemes (ISg 
RG, CSi RG, WH RG) and the Rayleigh-Schrodinger perturbation expansion (RSPE) restricting ourselves to the 
second-order perturbation expansion in powers of the bare quartic coupling qb- While the first-order contributions 
and the sign of the second-order ones are reproduced by all of these approximations there has been found a remarkable 
scheme-dependence of the weights of the second-order terms. 
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Appendix A: CSi RG scheme for the iV-dimensional anharmonic oscillator in real-time formalism 

1. Neumann-expansion 

For later use let us introduce the matrices 

Cw,w' — C(c,«),(c,«') = G u ',u, (Al) 



and 



B u ,-u = j [v(u)i,-ui) +v(u,-u})]q aiUl q a - Ul = -2 j ^~y(ul + w 2 ) + yQ g , Wl 9o - Ul , 

C w - W = -4^ [v{ui,-u) + v(-w 1 ,u)\q a , ul q a - ul =-% J (-^y^u + ^Y(uj 2 + oj 2 ) + q a , Ul q a _ Wl . 

(A2) 

The first two terms in the expansion (14) are quadratic in the field variable g 0)W and can be expressed in terms of the 
loop integrals I n ^ s defined in Eq. (17) in a straightforward manner, 

~ / G - = -y 7 ^ ( A3 ) 



2 y g ! w- | 4(o,w),(o,-w) — — y G^Bu-u + - j g^cuj 

= i[7Vy- (7V + 2)y]/ 2i0 / w 2 ga , Wl g a ,_ Wl 

+ 1 ([JVy - (TV + 2)Y] J 2 , 2 - {N + Q 2)9 J 2 ,o) y ?a >Wl ga,- Wl . (A4) 
The third term of the series (14) can be rewritten as 

GwA( c ,w),(d,u')Gu'A(d,-u'),(c,-u) 
= — / G 2 G u '[JVBui,i<j'-B- u - u ' + 2£?_ W; _ w 'C WiW / + ^( c ,w),(d,w')C'(c,-w),(d,-w')] 

^ Ju!,oj' 

= TX+T2 + T3, (A5) 
where all terms are quartic in the Fourier-transform of the field variable, 

T\ = — — f Gj.Gu'Buw'B 



2 

+2v(uj, J)v{ut ,uj 2 ) + v 2 (uj, oj')j q a ^ q a ,u> 2 qb,u: 3 qb,u 4 , (A6) 

— ill . — ill' 



12 - ~ fl ^ U ^«'^-u,,-u,'l> U)U 



J 

J LJ.LJ' ,UJl,LU2 ,W3,W4 V 
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T 3 



,w),(d,u/) C(c,-«),(d,-U)') 



— — 8/1 / G ! ^G ! w '(5 Wl +...+ W4j o5 Wl + W 3+ w + w ',o 

x [w(<x>i, cj)w(cj 2 , — w) + v(lji,u)v(uJ4, —oj') + v(uj 3 , oj')v(c02, — w) + w(w3, w')w(w4, —a/)] 



(A8) 



Evaluating 7\, T 2 , and T3 we keep only the terms corresponding to the interaction vertices included in the Ansatz (6) 
for the effective action in the truncated gradient expansion, ie. we keep self-interaction terms up to the quartic ones 
and quadratic gradient terms with wave-function renormalization truncated at the quadratic field-dependent term. 
By making use of the Dirac-deltas expressing energy conservation in the vertices, G u i can be changed to G u+a with 
a = ±(wi +W2). In order to reveal the terms of the gradient expansion explicitly, we have to expand the propagators 
with shifted frequency G u + a in the Taylor series (16) at a = 0. Furthermore, we also make use of the identity (9), 
in order to transform all gradient terms in one of those of the Ansatz (6). A rather lengthy but straightforward 
calculation yields finally, 

T i = ~ 2N J ^u 1 +u 2 +u 3 +u i fi^y 2 hA - 2-^yh,2 + 4 J^ / 3,0 



+ 



2 

2y 2 x 2 I 5 . 6 - 16^yx 2 I 5A + 32^x 2 h, 2 



~\xy 2 U A + 12^xyl 4 . 2 - &-^xI 4 ,o 



CO 2 



+ 



2y 2 x 2 I 5 . 6 - m^yx 2 I 5A + 32^z 2 J 5 , 2 



9 2 2 , 



24 2 ' 



-^xy 2 h A + 12-^yxl 4 , 2 ~ 8-^xI 4 ,o 
-y?h,2 + Y 2 h, 2 + 2^-yI 3 .o + 2^-YI 3 , 



24- 



24 



u 2 Qb,u 3 Qb,u 4 , 



(A9) 



To = -8 



+ 



-\yYh A + 2|?/ 3 , 2 - f 4 yh,2 + 4^/3,0 



-2x 2 yYI 5 . 6 - 8j^x 2 yI 5A + Wj^x 2 YI 5A + 32j^x 2 I^ 2 



+ ^xyYI iA + 0>i^xyh.2 - 12±xYh, 2 - >^_.rl 



+ \fha - h\yYh a + \y 2 H,2 + ^YI 3 .o - ^yh. 

9 2 2, 



+ 



-2x 2 yYh M - 8^.r 2 ///,.,, - 1 , ,- W - 32^,^., 



+\xyYI lA + ^xyh.2 - \2^xYU.2 - S^.rJ, () 



9_ 2 

24 X 1 ' 5 ' 4T ^24^ 
2 

24~ 2 



+ 4^/3,2 - 3-yr/ 3 , 2 + Y 2 ha + 2^^2/3,0 



LdiLU 2 



(A10) 



24 



and 



Jlo-i 



/ <Ll + ...+W4,0 

UJ\ ,...,U>4,U> 

Y\ gY T a 2 
T M + 12 ' 2 + 242 3,0 

+ (Vvj 5 , 6 + s^x 2 i 5A + 16^2^/5,2 



r 2 ff y o 2 



1 



52/ 



,5? 



+y 4,2 - -jV 4,2 - -j-h,2 - 24/3,0 + S^/S.O J^i 

( -y 2 ar 2 / 5 ,6 - 8^x 2 / 5 .4 ~ 16^^/5,2 



24 



24 2 ' 



Y 2 gY g 2 

+5— x/4,4 + 6—xU a + 4^x74,0 

3-2, 1^2 r , , 3£ r 1 

— -^y ^3,2 — 2 3 > 2 + ~4~ 3 > 2 + 24 — "24" l WlW2 



(All) 



2. Loop integrals 



In order to get the explicit forms of the RG equations (19)-(24) one has to perform the loop integrals (17), 

r0 ° du 



r°° du; r 



2irx n {lo - a) n {uj + b) n 



(A12) 



For the limit A — > 00 and for s — 2n < — 1, like in our case, this can be done analytically. Since Re a > 0, Im a > 0, 
we can close the path of integration with a half-circle of infinite radius on the upper half of the complex w-plane 
encircling the multiple pole at lo — a and apply the residuum theorem, 



1 

27 



(z-a)« (n-l)!da»- l/W ' /W (a + &)' 



Thus one finds 



4,o 



4,s 



i[s- (2n-3)] • (s - 1) 
(n - l)!2"x £ * i (r2-^ 2 ) n -T : 



— for n > 1, 



which yield the relations 



113 11 

4,2 = 2~^^ 1 - ' ^ 3 ' 2 = 4^^ 2,0 ' ^ 3 ' 4 = Sx 2 ^ 1 ' ' ^ 4 ' 2 ~ 6x^ 3 ' ' ^ 4 ' 4 = Sx 2 ^ 2 ' ' 



1 

8x~ 



1 



4,2- -4,0, 4,4-^- 



3,0, i 5,6 



64x 3 



'2,0 



with 



'2,0 



4xi (Q—i/j 2 )? 



, 4,0 



3i 



16x2 (Q— i\x 2 )\ 



, 4,0 



-hi 



32x^(^-^ 2 )2 



(A13) 



(A14) 



(A15) 



(A16) 



3. Perturbation expansion 



Here we outline the solution of the evolution equations (25)-(29) for the N = 1 dimensional oscillator in the second 
order of the perturbation expansion in the bare coupling gs after the analytic continuation — i/1 2 — > A to the Euclidean 



25 



space. Making use of the result of the IMA, one can write the perturbation expansion of the couplings as, 

1 = 5 + g B a + g 2 B l3 + . . . , 
Sfl 2 = g B p + g B a + 

x-1 = fflC + ■ • ■ , 

T = g%v, 

9 = 9B+9 2 B v (A17) 

with <50 2 = ft 2 — luq, and using x — 1 ~ O (T) <~ O (g%)- For the sake of simplicity, we shall set h = 1 in this section. 
In order to make explicit the dependence on the bare coupling g B in the right-hand sides of Eqs. (25)-(29), one has 
to expand the denominators of the loop integrals given by Eq. (A14), 

1 1 9bP 2 (\ , .-.I P 2 aa \ 

(f)2 + A)« K + A)° a K 2 + A)«+i + ' 9s ^2 ala+ \ Wo 2 + A)«+2 (wg + A)«+V' 

x a w l + ag|C + ..., (A18) 

where the relation 

(1 l x)a ~ l-ax+i a ( a + l) 2 ; 2 + ... (A19) 

has been used. The comparison of the coefficients of the corresponding powers of g B in both sides of the Eqs. (25)-(29) 
yields the following set of coupled ordinary first-order differential equations for the running couplings, 

= ikW< (A20) 

9x13 = 32( Wo 2 + A)l + 8(^ + A)l + 8( w g + A)5' (A22) 

^aC = ~ - ^— r , (A23) 
4(cj 2 + A) 5 

^ = 2^U 7 ' (A24) 

64(wg + A) 2 

Sap = - / 3 , (A25) 
8(a;g + A) 2 

dx<J = ~ 77 J 1 ~ — ^ rrr + ^ 7 ? ^ > ( A26 ) 
4(cj 2 + A)2 8(u; 2 + A)2 16( W 2 + A ) 5 
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The solutions belonging to the initial conditions 8b = — ^(w 2 + ^b) 2 , ctB = Pb = Cs = Vb = Pb = fa = Kb = at 
As — > oo and their IR limits A — > have been found analytically, 



'■/ 



1 



5 ' 



32( W 2 + A)f 32w§ 
1 1 

P ~ 4(w 2 + A)^ W 

3 3_ 

~8(o;g + A)7 

11 11 



256(w 2 + A) 2 256w£' 
1 1 



c 



32 {ujI + A) 32w 2 ' 

1 1 



3 ■ 128(w§ + A) 3 3 • 128wg' 



71 71 
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Then one obtains for the energy of the ground state 



Eq = 70 = -wo 



1 + «-i« 2 



for the frequency parameter 

the renormalized quartic coupling 



(A29) 

=o;g(l +4^- lie 2 ), (A30) 



9o=9b(1-60, (A31) 
and the couplings for the wavefunction renormalization 

2 



where we reestablished the powers of h. 



^o-l + |C 2 , fiT = 8^ 2 , (A32) 



4. Localization 



The localization of the ground state wavefunction can be characterized by its momenta. We can make an estimate 
of the dependence of the variance of the coordinate operator on the dimension N of the oscillator. Let us write for it 
with the help of the Lehmann-expansion of the propagator, 



N N 



<(% 2 |0) = ^(0| g 2 |0)= t lim + ^^(0|g Q (t)|n, a )(0,a|g a (0)|0) 

a— 1 a— 1 n, a 

= t Hmf;^e^- B «)*/ ft K0| go |n,a)| 2 

a— 1 n,a 
ii I n.a 



where H\n,a) — E n \n, a) for the exact eigenstates \n,a) and energy levels E n of the Hamiltonian, |0) and Eq stand 
for the ground state and its energy, a are additional quantum numbers for counting states belonging to the same 
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degenerate energy level. Operators without and with the time argument stand for those in the Schrodinger and the 
Heisenberg representations, respectively. Making use of the dominance of the pole of the integrand corresponding to 
the first excited state [27], one finds the estimate 

= ^E^E 2 ^™-^)^°I^K a )| 2 - (A34) 

a— 1 n,a 

Making use of the Thomas- Reiche-Kuhn sum rule [40], valid for any self-adjoint operator A, 

^(0\[[M,H],M}\0) =J2(En - E Q )\(0\M\n,a)\ 2 (A35) 

and the form of the renormalized Hamilton-operator H = ^xp 2 + V(q) where we neglected the less important higher- 
momentum terms in respect to the small values of y and Y, we get 

J2(En-E a )\(0\q a \n,a)\ 2 = \ (0| [[q a , H], q a ] |0) = |<0|[[« o ,^],« o ]|0) = \xt? (A36) 

n.a 

where there is no summation over the index a. With the help of this relation we obtain 

a=l 

Appendix B: WH RG flow of the Wilsonian action in Euclidean space 

1. Blocking via mode-by-mode integration 

The WH RG method [29] has been applied to simple quantum mechanical systems with success in the literature 
[27]. The simplicity of the + 1 dimensional model of the single quantized anharmonic oscillator enables one to 
perform the blocking by integrating out the high-frequency quantum fluctuations mode-by-mode. Let us prescribe 
periodic boundary conditions in time, q t= o — <Zt=T, £*=o = £*=t both for the coordinate q t and the mode £t integrated 
out in a single blocking step, where T stand for the time interval the action is integrated over and taken to infinity at 
the end. The Matsubara frequencies are given as u) n = n2w/T, the gliding cut-off is at ujn — N2ir/T. The periodic 
coordinate can be Fourier expanded as 

N A N a 

ft = y + ]>>n cos(o;„i) + ^ b n sin(w„i) = q + r\ t (Bl) 

n— 1 n— 1 

where ujn a = N\2n/T = A is the UV cut-off. Integrating out all the quantum fluctuations r\ t around the constant 
mode q=\a§ leads to the effective action So(q) given by 

e -\s {q) = J Vrje -\ S A[q+v\ 

/N A poo poo \ 

H da n db n \e-^+^ (B2) 

\ n —l^—oo J — oo / 

up to an irrelevant constant being independent of q and the potential (sec [42]). The blocking relation for integrating 
out the single mode with frequency ojn, ie. the mode £,t — <^n cos(uJNt) + &iv s'm(ujNt) is given by 

/roo poo 
V ^ e -is N [ q+S ] / daN dbNe - iSN [q +i ] (B3) 
J — oo J —oo 
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in the LPA. Here the LPA enables one to replace the time-dependent background by the constant mode q. The kinetic 
term of the action gives 



(B4) 



dt\{qo+it) 2 = ^ 2 N [a 2 N + b%}. 
With the notation a a= \ — ajy, a CT =-i = frjv ; we can write for the blocked potential 

U N [q + Z\ = fu N {q + i t )=TU N {q)+ l - ]T a a a a ,U (2) ^, (?) + O (f) , 



where 



so that one gets 



= T^w 2) (9)[^,i^M + k,-i$x',-i], 



+ = TC/ Ar (g) + r i c/^ 2) (g)[ a 2 r + 62 r ]+o(e 3 ). 
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The non-Gaussian higher-order terms are suppressed with some powers of \/T in the path- integral and can be 
neglected. Then the Gaussian path integral can be performed, 

da N / db N e~i T ^- +u - 

-oo */ —oo 



~ e -*™«(<7) 



8irh 



T[u;% + U%\q)} 

yielding the finite difference equation 

U N (q) - UN-i(g) = Uu, N {q- Uu N -Au(q) 

2tt/T Aw 2tt/T' 

for the blocked potential at the scale con- In the limit T — > oo, ic. Aw — > one finds the WH RG equation 



/,/T ; m[<4 + £/£%)] 
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(BIO) 



for the blocked potential U u , where the discrete gliding frequency cut-off ojn has been changed to the continuous one, 
denoted by cj. 



2. Approximation schemes 

a. IMA 



In the IMA one has to replace the blocked potential in the right-hand side by the bare one. Then one can perform 
the integration over the scale ui in a straightforward manner, 



UT A (q) = U A (q) + Ha f dk\n[k 2 + a 2 ], a 2 = U {2 \q). 



(Bll) 



The bare potential is convex at least in the neighbourhood of its minimum (minima) at some qo- Then for q sufficiently 
close to qo, one can make use of the integral 



,A 

/ dkln[k 2 +a 2 

J UJ 



kln[k 2 + a 2 ] - 2k + 2a arc tan 
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(B12) 
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and find for the effective potential 

U™ A (q) = U A {q) + tmy/u^{q). (B13) 

This expression is valid in the neighbourhood of the minimum qo of the bare potential, and the UV divergent piece 
is removed by the renormalization condition that the effective potential of a free particle (J7a(?) =0, a = 0) were 
identically vanishing. 

The minimum at q of the bare potential would be displaced to q$ MA = qo + Sq with 

2pf (»)l s/2 

so that for a potential with the symmetry q — > —q and its minimum at qo = the minimum is not displaced in the 
one-loop order. Then one finds for the ground state energy 



El MA = hcxxyjuj?\0) (B15) 



in the IMA, ie. in the one-loop approximation. For the free motion of a particle the bare potential vanishes and it 
remains vanishing after the blocking as well. For the case of a linear harmonic oscillator this reproduces the correct 
ground state energy i?J MA ' lho ' = ^huoo- For an arbitrary polynomial bare potential U\(q) = | w o9 2 + ^i9Bq 4 + ... we 
get 

U™ A > aho -(?) = t/ A (,-) + ^ (l + ^ + ...) = i^ + ^(a; 2 + g) 9 - 2 + ^W + ... (B16) 

This means that the energy of the ground state of the anharmonic oscillator remains the same as that of the harmonic 
oscillator, but the excitation energy of the first excited state is increased with increasing coupling g B , 



9Bft f. /, . 9b^ 
- — wto 1+ — -3 

in the IMA. The same result has been obtained in the CSi RG scheme in the IMA. 



a/:- - = f g « fi^o^i + g£ 4 o (ni) ) (bit) 



"(fll)). 



6. Perturbation expansion in powers of the bare coupling gs 
Let us introduce the blocked interaction potential V u (q) via 

U u (Q)=u u + ^lq 2 + V u {q), (B18) 
for which the WH-RG equation (BIO) can be rewritten as 

00 —(—\) n 

d„u u + 8Mq) = -haHu: 2 + u 2 + V^{q)] = -ha\n{uo 2 + u 2 ) -haj^ ^^[G ^{q)\ n . (B19) 

n 

n=l 

In the last equation we performed the perturbation expansion of the right hand side into the interaction, Go denotes 
the propagator of the harmonic oscillator in the ground state. Further on let us Taylor-expand the blocked interaction 
potential into the powers of q and truncate it at the quartic terms, 

V u {q) = gBa^ + g 2 B /3^ + ^(g B p^ + g 2 B a^)q 2 + -^(g B + g^^q 4 + . . . , (B20) 

and Taylor-expand each running coupling into the powers of the bare coupling g B and truncate their perturbation 
serii at the second order terms. Inserting this ansatz into Eq. (B19) and keeping all the terms up to the order q 4 
and g 2 B , one finds the following set of coupled ordinary differential equations for the RG flow of the various couplings 



30 



(' = d/dui), 



O-u. 

k 
Pu 



ftaln(w 2 + uJq), 

Pu 



—ha 
—ha 
—ha 
—ha 
3fta 



to 2 


1 




<*U 


9 
W 2 


+ 0J 2 




1 


w 2 


+ < 




Vu + 


W 2 






1 



lha-A 



,2^2 ' 



ha- 



,2\2 ' 



(w 2 +W 2 ) 2 



(B21) 



with the initial conditions as — Pb = Pb = ob = = at the UV cut-off w = A. Here the equations for 
Uu, Pu, Vu-, &u can be integrated straightforwardly, 



,A 

Uu = u\ + ha dk ln(fc 2 + w 2 ) 

J u 



A— 7-oc. u— 7-0 TTT7 - ,. 1, 

— ■> UV div. terms + -ftw , 



(B22) 



Pu; = 



2tt 



fc 2 + w 2 27rw 



arc tan 



w 



A^ oo, u ^Ct h 

4w ' 



(B23) 



3ft r A 



dk 



(fc 2 +w 2 ) 2 



3ft 



k 1 arctan k 

2w 2 (fc 2 +Wg) 2uj^ w 



n A 



A^oo, u^O 



5' (B24) 



ft 2 a 2 f A arc tan — — arc tan — 



w 
fi 2 or 



/ 



fc 2 + w 2 



w 

ft 2 a 2 r 



1 A / fc 

— arc tan — arc tan — 

w w V w o 



1 

2wo 



A „ J k Y k 

dk 2 1 arc tan — arc tan — 

wo/ w 



w 



1 A , 

— arc tan — arc tan 

w w V w o 



*Y-J-f 

2w V 



arc tan — 

w 



A^oo, w->0 h 2 



32w§' 



(B25) 



The rest of the equations cannot be integrated straightforwardly, but their solutions can be approximated by inserting 
the IR limits of the couplings p u , into their right hand sides, 
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Finally one finds 
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= u +g B u +g 2 B l3 {) = -hw {) {l + Z-l%f + ...), 
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ft 2 = w 2 + . 9B/ oo + ffl^o = w 2 (l + 4£ - 32£ 2 + ...), 
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9b + 9 B vn = 9b{1 - 6£ + . . .) 
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in terms of the dimensionless parameter £ = gsft/(16wg). One should mention that the gradient expansion cannot be 
improved going beyond the LPA in the framework of the WH RG scheme. 



